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Abstract 

A supersymmetric way of imposing the constraint of no double occupancy 

in models with strong on-site Coulomb repulsion is presented in this paper. 

In this formulation the physical operators in the constrainted Hilbert space 

are invariant under local unitary transformations mixing boson and fermion 

representations. As an illustration the formulation is applied to the t — J 

model. The model is studied in the mean-field level in the J = limit 

where we show how both the slave-boson and slave-fermion formulations are 

included naturally in the present approach and how further results beyond 

both approaches are obtained. 
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The t — J model has became a focus in the study of strongly correlated metals and 
High-T c superconductors since it was proposed in late eighties [|I|]. Because of lack of small 
parameters for expansion, analytical understandings of the model were largely depending 
on mean-field theories which treat the constraint of no double occupancy only on average. 
So far, the most successful mean-field approaches to the t — J model seems to be based 
on either the slave-fermion mean-field theory(SFMFT) [|J which is successful at very small 
doping when antiferromagnetic correlation is important, and the slave-boson mean- field 
theory(SBMFT) which is successful at larger value of doping when the system becomes 
superconducting. The only difference between the two approaches is that two different 
representations of spin and electron operators are used to impose the constraint of no double 
occupancy. More recently, the focus in the study of High-T c superconductors has turned to 
the underdoped and spin-glass regimes where it is believed that the subtle interplay between 
antiferromagnetism and superconductivity determines the properties of this crossover region. 
In particular, the importance of SU(2) symmetry in the underdoped regime of the t— J model 
has been pointed out Alternatively, it was also suggested that an 5*0(5) symmetry 

may play an important role in determining the competition between antiferromagnetism and 
superconductivity in the high-T c cuprates ||. 

To understand the complicated behaviour in this regime of the t — J model, it seems that 
a unified approach which incorporate both the advantages of the slave-fermion mean-field 
theory and the slave-boson mean-field theory is essential. In this paper, we shall show that 
it is possible in general to formulate models with constraint of no double occupancy in a 
way which incorporates the advantages of both slave-fermion and slave-boson representa- 
tions. In this new formulation the physical operators are super symmetric and are invariant 
under unitary transformations mixing fermion and boson representations. The formulation 
suggests that supersymmetry exists naturally in strongly-correlated systems where on-site 
Coulomb repulsions are strong. In the following, we shall use the t — J model as an example 
to illustrate our approach. To begin with, we first consider the Hilbert space of a lattice 
model with constraint of no double occupancy imposed. 
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The constraint of no double occupancy implies that there are three possible states on 
any single lattice site in the model. The site can be either empty (hole state), or can be 
occupied by either an up- or down- spin electron. In the slave-boson approach, the hole state 
is represented as a boson, whereas spins are represented as fermions 0. It is also equally 
valid to represent spins as bosons, and holes as slave fermions, as in the slave fermion 
treatment ||. In our formulation we shall consider an enlarged Hilbert space where both 
possibility of representing hole and spins coexist as different states of the system, i.e. there 
are now six possible states per site in the Hilbert space, represented by, 

>= c+ |0 > , \hP >= 6+|0 >, (la) 

and 

|of >= z ia \o > , \h!p>=f?\o>, (lb) 

where a =T, !, cf a and Z ia are fermionic and bosonic spin creation operators, respectively 
and |0 > is the vacuum state. Similarly, bf and are bosonic and fermionic hole creation 
operators, respectively. Notice that we have seperated the states into "slave-boson" (la) 
and " slave-fermion" (lb) groups in Eq. (1). For a system of iV-lattice sites, both groups 
of states are allowed at all sites in our formulation and the total Hilbert space is thus 2 N 
times larger than the Hilbert space of the original model. The essence of our approach is to 
construct a Hamiltonian which is equivalent to the original model in all these 2 N groups of 
states, and consequently our system with enlarged Hilbert space is equivalent to 2 N replicas 
of the original model. To see how this Hamiltonian can be constructed for the t — J model 
we first consider spin operators. 

We consider the spin operator s, at site i, 

Si = l? ) +Si 6> ) (2a) 



where 
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and a is the usual Pauli matrix. Notice that sU> and s^ b > are the spin operators in usual slave- 
boson and slave-fermion representations, respectively. It is obvious that ifj is itself a spin 
operator since it is the sum of two spin operators. The matrix elements < <j\ \~Si\o i > 
where a, /3 = f,b and a, a' =|, I can be computed easily where it is easy to see that 
< af^KP >=< vi b) WiWf > and gives the usual spin operator matrix elements between 
spin-1/2 states whereas all other matrix elements with a 7^ (3 are equal to zero. Thus the 
states \a\^ > and \af^ > together form two identical replicas of spin-1/2 states on site i 
with our definition of spin operator (0). It can then be shown by direct evaluation that the 
Hamiltonian 

Hj = J Si.sj, 

<i,j> 

represents 2 N identical copies of Heisenberg interaction in our system of enlarged Hilbert 
space. 

The electron annihilation and creation operators in our system can be defined as ipi a = 
hf^ia, and ij)f a = ^hi, respectively, where 







( z- \ 















are doublets of hole and spin operators carrying fermion and boson statistics. It is straight- 
forward to show that the electron operators defined this way satisfies the usual electron 
commutation relations in the Hilbert space spanned by Eq. (1). The kinetic energy term 
H t of the t — J model can be constructed by requiring that the hopping matrix elements 
< erf, h 1 ^ \H t \hj,a^ >= —tij8 aa i8pj3', where a , a' , (3 , f3' = b, / and \h",a^ > represents a 
state with a hole belonging to group a on site j and a spin o belonging to group (3 on site i. 
Notice that the group indices a, /3's are "conserved" in constructing H t . It is straightforward 
to show that 

H t = ~t ( C t b i b t C ™ + C^bjffZia + Zjjjbfcie - ZjafjffZia + C.C.) , (4) 

<i,i>,<7 
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where we have considered a Hamiltonian with nearest neighbor hopping only. The four 
different terms in H t give the matrix elements of the hopping term between the four possible 
combination of groups of states at sites < i,j >. With this it is easy to verify that our 
system with Hamiltonian H = Hj + H t is equivalent to 2^ copies of the usual t — J model. 
Notice that the hopping term cannot be simply represented as H t = — tXXVwrVv + c - c -) 
because of the sign difference in the " slave-fermion" term Zj a ]jjfZi c . This sign difference 
is well known in studies of slave-fermion mean-field theory ||. 

The invariance of our system under change of group of states defined in Eq. (1) at any 
lattice site can be expressed in the language of supersymmetry. To see that first we note 
that the spin operator (0) can be written using the £ CT fields as 



where f£f i(7 ' = Z ia Z ia < + 

c ia c ia- It is now easy to see that the spin and electron operators 
are invariant under the super-unitary transformation 

Cia -> U& a , hi -> Uihi (6) 

where U^s are local 2x2 unitary super-matrices mixing the fermion and boson representa- 
tions of spins and holes. The generators of the super-unitary transformations are superspin 
operators 

st = \ ((E c i c - + W - (E + ftfJj 

St = E C taZia —b^fii 
Si — E Zia c ia ~ /i + ^l> 



(7) 



where we also define the (magnitude) 2 of the superspin as S 2 = S Z S Z + ^(S + S~ + S~S + ). 
Using the fact that the allowed states in our Hilbert space can only be singly ocupied by 
either spin or hole, it is easy to show that 



(8) 



where we have used the requirement that J2a(Zi a Z ia + cf a Ci a ) + ft fi + bfbi = 1 in our 
Hilbert space (1) in writing down the last equality. Notice that the constraint of no double 
occupancy is equivalent to the condition that the magnitudes of the superspins are fixed 
(= 1/2) on all lattice sites! In terms of the superspin operators, the slave-boson and slave- 
fermion representations are equivalent to fixing the direction of superspins to be pointing up 
and down, respectively and supersymmetry expresses the fact that the physical observables 
in the system are in fact, invariant under local (but time-independent) rotations in the 
superspin space. Notice that because of the "minus" sign in the " slave- fermion" hopping 
term Zj a fjftZ ia , supersymmetry is straightly speaking, broken by the t— term in the t — J 
model. However, our formulation of t — J model in the enlarged Hilbert space (1) does not 
require supersymmetry in the Hamiltonian and is still valid. 

The Lagrangian of the t — J model in the enlarged Hilbert space is 



and a Path-integral formulation of the problem can be written down as usual. In particular, 
the partition function of the present model is 2 N times the partition function of the original 
t — J model. Notice that the constraint of no double occupancy in our enlarged Hilbert 
space is imposed by an Lagrange multiplier term as usual. 

It is obvious that the same approach can be applied to other models with constraint of 
no double occupancy, as long as a proper Hamiltonian can be constructed in the enlarged 
Hilbert space. For example, we find that the Hamiltonian for the Infinite-C/ Anderson model 
in the supersymmetric representation is, 

Hand = e k a t a ka + £ »E + V ^(a^Ur + ^£a«r), (10) 

where Q>{o> + )% c is the annihilation(creation) operator for conduction electrons with momen- 
tum k and spin a and £(£ + )i CT are annihilation(creation) operator for localized spins on site % 




j 



(9) 
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where constraint of no double occupancy is imposed, eg, e Q and V have their usual meaning 
in Anderson model. The spin and electron operators in the localized orbitals are represented 
by £«cr and operators which are supersymmetric in our formulation and the constraint 
of no double occupancy can be imposed by Langrange multiplier fields as in the t — J 
model. Notice that the Infinite-t/ Anderson model is, in fact, completely supersymmetric, 
in constrast to the t — J model where supersymmetry is broken by t— term. 

To see the advantage of present formulation we shall consider in the following the t — J 
model in the J = limit, and shall study the model at two dimension in the mean-field 
level. The model is equivalent to the U — oo Hubbard model which is a model of interests 
in itself |7|f|. 

A mean-field theory in the supersymmetric t — (J = 0)-model can be obtained by making 
the following decouplings, 



h -> -t {< 4^ > b i b t+ < h M > c t* c i°- < c l* c i* >< h M > + c - c - ( n ) 

<i,j>,a 

+ < c^Zfc > b 3 ft + cj a Z ia < b,jf + > - < c^Ziv >< bjf + > +c.c. 
+fjbf < Z ja c ia > + < fjbf > Z ja c ia - < fjbf >< Z ja c ia > +c.c. 
— < ZjfjZicr > fjfi~ — Zj a Z i(J < fjfi~ > + < Zj a Z i(J >< fjfi~ > +c.c. ^ 
+A 52(Z ia Z ia + ctc la ) + (A - ai) + Kh) 

i,cr i 

where < A > is the expectation value of operator A evaluated with the mean-field Hamil- 
tonian. The mean-field parameters A and \x are chosen so that < ft fi + bfbi >= 5 (hole 
concentration) and J2a < Z i(T Zi a + c[ a Ci a >= 1 — 5. We shall be interested at translation- 
ally invariant solutions where the mean-field parameters are independent of < i,j > in the 
following. 

Using the fact that < C >= for Grassman variables C's, we find that the second 



and third lines in the mean-field Hamiltonian (|TT| ) are equal to zero. However, both slave- 
boson and slave-fermion type terms are still present in the mean-field Hamiltonian (first 
and fourth lines in flTTp). The relative weight of the two kinds of mean-field terms are 
determined by the occupation numbers of spins (and holes) with fermion (boson) and boson 



(fermion) statistics. These numbers are in term, determined by minimizing the free energy 
of the system with the constraint that the total number of spins and holes are fixed to be 
1 — 5 and 5, respectively 

Solving the mean-field equations we find that at low temperatures there are two solutions 
corresponds to local minima in the mean-field free energy for any given hole concentration 
5. The two solutions have either < c + c >, < b + b >^ and < ZZ >, < f + f >= or the 
other way around and corresponds to the usual slave-boson and slave-fermion solutions. New 
solutions where both slave-fermion and slave-boson mean-field parameters are nonzero are 
also present. However, these solutions correspond to saddle points in free energy landscape 
and are unstable. Comparing the free energies of the two stable solutions we find that at 
zero temperature and at low doping 5 < 0.33 the slave-fermion-like solution has lower energy 
whereas for S > 0.33 slave-boson-like solution has lower energy. The slave-fermion solution 
where the spinons condensed with < Z^ >=< Z^ > corresponds to ferromagnetic state with 
spin pointing in x— direction whereas the slave-boson solutions correspond to paramagnetic 
state. The mean-field theory predicts a first-order transition at zero temperature from 
ferromagnetic to paramagnetic state as concentration of hole increases across ~ 0.33. Notice 
that numerical J7[ and analytical || works have established that the Nagaoka (ferromagnetic) 
state is unstable in the infinite-U Hubbard model for any finite concentration of doping 5 and 
our mean- field phase diagram is incorrect. Nevertheless, our mean- field theory does produce 
the qualitative features that the spins (holes) excitations are bosonic (fermionic) like in the 
ferromagnetic state of the model, and are fermionic (bosonic) like in the paramagnetic state 

It has to be emphasized that although the mean-field solutions are either slave-boson 
or slave-fermion like, spins and holes excitations with both statistics are present in the 
supersymmetric mean-field theory. For example, in the slave-boson-like solution where < 
ZZ >, < / + / >= 0, bosonic spin and fermionic hole excitations still appear as dispersionless 
high energy excitations with energies A and A — /i, respectively in the mean-field theory. 
Similar situation occurs also in the slave-fermion-like solution. The simultaneous appearance 
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of spin and hole excitations with both statistics in the supersymmetric mean-field theory 
implies that at finite temperature, the physical properties of the theory is quite different 
from usual slave-boson or slave-fermion mean-field theories. For example, extra incoherent 
parts in the one-electron Green's function which does not exist in usual slave-boson or 
slave-fermion theories will appear in the supersymmetric theory at nonzero temperature. 

At high temperature T >> 5(1 — 5)t the only stable mean-field solution which exists is 
with mean-field parameters < ZZ >, < / + / >, < c + c >, < b + b > all equal to zero, indicating 
that the high-temperature phase is completely incoherent. The transition from coherent low 
temperature states to incoherent high temperature state is a first-order phase transition in 
mean-field theory. However, similar results were also obtained in conventional slave-fermion 
or slave-boson mean-field theories where it is believed that the prediction of first-order phase 
transition from low-temperature to high-temperature phases is a defect of mean-field theory 
and should be replaced by a smooth crossover when temperature increases. The mean-field 
phase diagram is shown in figure 1, where the three stable phases are seperated by lines 
indicating first-order phase transitions in mean-field theory. 

Summarizing, by extending the size of Hilbert space, we present in this paper a new way 
of formulating t — J model where supersymmetry is inherent in the physical observables. 
The consequence of this hidden supersymmetry is exploited in the t — J model where a 
new, supersymmetric mean-field theory of the model is introduced and is studied in the 
J — > limit. The new mean-field theory unifies the slave-boson and slave-fermion mean- 
field treatments of t — J model by including them as subsets of possible solutions of the new 
mean-field theory. More generally, new phases where both slave-boson and slave-fermion 
type mean-field parameters are nonzero may also exist. 

Perhaps the more important message in this paper is the demonstration of existence 
of supersymmetry in general strongly correlated systems where the low energy physics is 
described by effective Hamiltonians with constraint of no double occupancy like the t — J 
or infinite-?/ Anderson model. The consequence of this supersymmetry is never exploited in 
studies of strongly-correlated systems. Note that supersymmetry is broken by the t— term 
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in the t — J model and it's consequences are not fully exploited in our paper presenting 
only the mean- field treatment of t — ( J = 0) model. Further studies of the t — J model is 
underway and the results will be presented in future papers. 

We acknowledges the support of HKRGC through grant no. HKUST6143/97P. 
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FIGURES 

FIG. 1. phase diagram of the tt J = 0) model in supersymmetric mean- field theory: (HT)-high 
temperature phase, (SB)-slave-boson-like phase, (SF)- slave-fermion-like phase 
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Phase Diagram 
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